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Abstract 

We start from a low-energy effective field theory for interacting fermions on the lattice and 
expand in the hopping parameter to derive the nearest-neighbor interactions for a lattice gas 
model. In this model the renormalization of couplings for different lattice spacings is inherited 
from the effective field theory, systematic errors can be estimated a priori, and the breakdown 
of the lattice gas model description at low temperatures can be understood quantitatively. We 
apply the lattice gas method to neutron matter and compare with results from a recent quantum 
simulation. 
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I. INTRODUCTION 



Lattice gas models play an important role in many areas of physics from fluid mechanics 
[l|, 3, ^ to quantum computing 0, ^, alloy mixing El to nuclear physics. In nuclear 
physics phenomenological lattice gas models have been used to model multifragmentation in 
heavy ion collisions and the thermodynamics of symmetric and asymmetric nuclear matter 

at nonzero temperature 0, y, H H Q, Q, Q Q, Q- 

In this study we attempt to 

broaden the scope of the lattice gas approach as applied to interacting fermions at nonzero 
temperature. We build a connection between lattice gas models and low-energy effective 
field theory on the lattice. While our main interest concerns interacting nucleons, our 
approach to lattice gas models should apply to systems such as trapped Fermi gases near a 



Feshbach resonance |l7l . 
are applicable. 
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22| |. where similar effective field theory descriptions 



For most applications of lattice gas models in nuclear physics, the coefficients of inter- 
actions are treated as adjustable parameters tuned to make the model realistic. Also the 
lattice spacing is often chosen so that the completely filled lattice corresponds with normal 
nuclear matter density, ~ 0.17 fm~'^. While one cannot argue with the successes of this 
phenomenological approach, there remain some fundamental questions. How do we know 
which interactions are needed to describe the important low-energy physics? How can we 
determine the coefficients of the interactions directly from binding energies and/or few-body 
scattering data? How can we do simulations at different lattice spacings while keeping the 
low-energy physics the same? 

In the full quantum theory these questions are answered by effective field theory. In the 
ow-energy limit, power counting schemes organize the interactions in order of importance 
26l |. Over the last few years effective field theory methods have been used 
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There has also been 



to study two and three- nucleon systems at low energy [27, 
recent prog ress in applying effective field theory to many-body nuclear lattice simulations 
0, 0, In this approach, operator coefficients in the effective lattice Lagrangian are 
matched to few-body scattering data, and the renormalization group describes how the 
operator coefficients depend on the lattice spacing. The resulting lattice action can then 
be simulated using standard lattice Monte Carlo methods to produce many-body results. 
In this study we attempt to bring lattice gas models into the framework of effective field 
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theory. We make this connection by means of a spatial hopping parameter expansion. 
Starting from any effective field theory of interacting fermions, we show how to construct 
the rules for a corresponding lattice gas model. We discuss the convergence of the hopping 
parameter expansion as well as the uses and limits of the lattice gas approach. In particular 
we clarify why it can describe a "classical" phase transition but not a truly "quantum" phase 
transition. As an example we consider low-energy neutron matter and compare with recent 
quantum simulation results js^ for the energy per neutron as a function of density. 



II. HOPPING EXPANSION 



Let us consider n f species of fermions and any operator V built up from the annihilation 
and creation operators. Let \9) denote the eigenstates of V, 



v{a\a) \e) = v{e)\e). 



(1) 



Next we consider a tensor product space that represents the fermions at two different loca- 
tions. For the moment we single out one particular species, j, and define the Hamiltonian 



H- = V(a\ a) (g) 1 + 1 (g) V{a^, a) - e 



(g) cij + cij (8 



Let us define the matrix element 



4^''^(/5) = K^i|®(^2|]exp 



If we expand in f3e we find 



z'^''-i/3) = exp[-/3iVie^) + V{e,))]. l + /5V/f'^^(/?) +0ii(3er' 



(2) 



(3) 



(4) 



where 

6', 9" 

and 



|(^'|a,|^i)|2|(^"|at|^^)|2 

+\{e'\a]\e,)\'\{e"\a^\e2)\^ 



■F [f3{v{e") + v{e") - v{e,) - v{e2))] (5) 



Fix) 



1 + x 



x^ 



m 



1 



(6) 



We now generalize this result to a general three-dimensional periodic lattice Hamiltonian 
with nearest-neighbor hopping and single-site interactions, 



In the following we use dimensionless lattice parameters. If a is the lattice spacing then the 
dimensionless mass parameter and inverse temperature are 



m = mphysa, 
1 



(8) 
(9) 

(10) 



Let |G) be a configuration of fermion states at each lattice site, 

\&) = (Si\0{n)), 

n 

Let us define 

^(/3,e) = (0|e-^^|e). (11) 

After applying the hopping corrections in ^ for each lattice site, fermion species, and 
dimension, we find 

"2 



z{P, 9) = exp 
+ 0{{ 



-/?J]n^(n)) 



n n n 

_ H j=l,nf 1=1,2,3 



1 + 



P \ r.9{n+l),e{n) 

2^1 



i/3) 



2m 



(12) 



We have assumed that the lattice is longer than three sites in each dimension. If the lattice 
were three sites long in some dimension then there would be terms at 0{{-^)^) which wind 
around the lattice. 

We can introduce a chemical potential by adding —fiN to H, where N is the total fermion 
number operator. In order to compute 



we can use the same expression in (jl2p if we redefine 

V ^ V — ^ Ojdj. 



(13) 



(14) 



By summing over all configurations G we now have an approximation to the grand canonical 
partition function. 



ZG = Tr [e-^(^--^)] =5^z(/5,/.,e) 

e 



^exp 







n n n 

n j=l,nf 1=1,2,3 



(15) 
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We note the similarity to a 2^-f-state Ising model. The only differences are the non- 
exponential factors and the temperature dependence in the interactions. One could write 
everything in exponential form using 



1 + 



1_ 
2m 



f 

J 3 



(n+l),d{n) 



{j3) ^ exp 



2m 



^{n+l),d(n) 



(16) 



For weakly coupled systems either form will do. However we find that for very strongly- 
coupled systems the exponentiated form produces larger 0((^)^) errors than the original 
expression. 



III. CONVERGENCE AND LONG-RANGE ORDER 

The spatial hopping parameter expansion can be extended to higher orders. At order 
0((^)") one must consider all n-step paths which are connected and form closed loops. 
On an lattice where L is even, all closed paths must have an even number of steps. In 
that case only even powers of ^ are nonzero. When L is odd, there are winding paths that 
give nonzero contributions for odd powers greater than or equal to L. A similar expansion 



was used to derive the zone determinant expansion in j33|, where the parameter ^ was 
identified as a localization length in lattice units for a given fermion. 

As the temperature increases the hopping parameter expansion converges more quickly. 
However if the temperature is too high then the relevant physics may be at momenta too 
high for the chosen lattice spacing. The momentum cutoff scale on the lattice is A^ = vra"^. 
The temperature must therefore lie well below the kinetic energy associated with this cutoff 
scale, 

Tphys « ^r^. (17) 



In order to have a sensible effective theory we need 

1 (3 
— ^ 



« (18) 



If we combine this with the condition for convergence of the spatial hopping parameter 
expansion we get 

1 ^ /5 _ 1 ^ ^ ^^^^ 

7r^ 2t71 2mip}iygTp}iygOp' 
Let us consider as an example T 20 MeV, roughly the temperature for the liquid-gas 

transition in symmetric nuclear matter. For a well-defined lattice gas model based on a 



hopping parameter expansion, (fT^ tells us that the lattice spacing must lie in the range 
from about 1 fm to 3 fm. This is enough to probe a wide range of densities, including 
the saturation density ~ 0.17 fm~^. Therefore it seems possible to describe the phase 
transition at several different lattice spacings. 

While a lattice gas model may describe long-range particle density ordering in a liquid-gas 
transition, it cannot describe long-range order associated with "off- diagonal" operators. By 
off-diagonal operator we mean operators which don't commute with the single-site operator 
V . In our lattice gas model formalism these operators are quite different from diagonal 
operators, such as the particle density operator, which commute with V . We can compute 
the correlation functions of diagonal operators simply by computing the eigenvalues associ- 
ated with each configuration state |0). But in order to compute the correlation functions 
of off-diagonal operators, we need to consider entirely new hopping paths connecting one 
operator to another. It is clear that any long-range correlations would have to be built by 
hand from arbitrarily long paths in our hopping parameter expansion. 

Therefore we expect the lattice gas model approach to be ineffective for any truly "quan- 
tum" phase transition. Long-range order in a quantum phase transition becomes possible 
only when the quantum wavefunctions of individual particles overlap. Hence the localization 
length ^ must be greater than the interparticle spacing in lattice units, and so therefore 
^ > 1 

2m ~ 

IV. APPLICATION TO NEUTRON MATTER 

We now apply our hopping parameter expansion to an effective field theory for dilute 
neutron matter on the lattice. We focus on low energies and densities where the relevant 
momenta are smaller than the pion mass, and we use an effective field theory with only 
neutrons. We work with the lowest order effective Lagrangian which contains a neutron 
contact interaction that is adjusted to produce the physical ^Sq scattering length. Our 
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lattice Hamiltonian with chemical potential included has the form 

n 

1 

n j=T,i«=l,2,3 

We will use the labels Q — 0,t,t,ti to represent the various zero, one, and two neutron 
states on a single site. In Table 1 we list V{Q., //) for the various neutron states. 

Table 1: V{e,^) 



~ ^ $Z $Z [«] (^) «j + + a] (n) a^- (n - I) 



(20) 
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In Table 2 we hst for the various neutron states on nearest-neighbor sites, and in 
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fi2, 



Table 3 we list 



Table 2: 








t 


i 


n 








1 

2 







T 


1 

2 










i 





F(/3C) 





1 

2 


Ti 







1 

2 






Table 3: 
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V. RESULTS 



We have run lattice gas model simulations for both free and interacting neutron matter. 
The value for the interaction coefficient C is set by comparing with experimental data from 
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nucleon-nucleon scattering. We sum all nucleon-nucleon scattering bubble diagrams on the 



lattice, locate the pole in the scattering amphtude, and compare wit 



relating scattering lengths and energy levels in a finite periodic box 
are shown in Table 4. 



iischer's formula 
The results 



1 u uau ii 

34, 3- 



Table 4: Interaction coefficient C for different lattice spacings 



a-i (MeV) 


C (MeV-2) 


50 


-8.01 X 10-5 


60 


-6.73 X 10-5 


70 


-5.81 X 10-5 


80 


-5.10 X 10-5 



We compute the energy per neutron, E/A, as a function of neutron density. The total 
number of neutrons. A, and average energy, E, are computed using 

1 d 



A 



d 



InZ 



(m — fi)A. 



(21) 
(22) 



The results for T = 8 MeV are shown in Fig. ^ In Fig. |21 we show a similar plot 
from a full quantum lattice simulation The lattice volumes for our lattice gas models 

. In both plots 



are chosen to be the same as that for the corresponding simulations in 
we use the abbreviation "fc" for free continuum results, "f ' for free lattice results, and "s" 
for lattice simulation results. In Fi g. El results for bubble chain diagrams calculations are 
also included and labelled with "b" [32]. In addition to these abbreviations, we also use the 
shorthand labels shown in Table 5 for various combinations of spatial and temporal lattice 
spacings. 

Table 5: Shorthand labels for various lattice spacings 



a-\MeV) 




Label 


50 


24 





60 


32 


1 


60 


48 


2 


70 


64 


3 


80 


72 


4 



For the lattice gas model however the temporal lattice spacing, at, is set to zero. 
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FIG. 1: Lattice gas model results for energy per neutron versus density at T = 8 MeV. 
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FIG. 2: Full quantum effective field theory results for energy per neutron versus density at T 
MeV. 



In Fig. 01 we show results at T = 4 MeV for the lattice gas model, and in Fig. ^we show 
the full quantum simulation at T = 4 MeV. For each of the temperatures and spatial lattice 
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FIG. 3: Lattice gas model results for energy per neutron versus density at T = 4 MeV. 
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FIG. 4: Full quantum effective field theory results for energy per neutron versus density at T = 4 
MeV. 



spacings we show the corresponding spatial hopping parameter in Table 6. 

Table 6: Hopping parameters for various lattice spacings 



a-i(MeV) 


50 


60 


70 


80 


for T = 8 MeV 


0.17 


0.24 


0.33 


0.43 


/- for T = 4 MeV 

2m 


0.33 


0.48 


0.65 


0.85 
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As discussed in the previous section, the lattice gas model cannot describe the superfluid 
transition in neutron matter since this requires long-range ordering associated with neutron 
pairing. There may be some indication of this already in the T = 4 MeV data shown in 
Fig. 121 The points at higher density show considerable deviation from the full quantum 
simulation results in Fig. 0J Nevertheless we see that the lattice gas results agree quite well 
with the quantum simulations when the spatial hopping parameter is less than about 0.4. 
This is a bit surprising considering that the computational cost for the lattice gas model 
simulation is several hundred times less than the quantum simulation. 

VI. SUMMARY 

Starting from a low-energy effective field theory for interacting fermions on the lattice, we 
derive the nearest-neighbor interactions for a lattice gas model by expanding in the spatial 
hopping parameter. Unlike most phenomenological approaches, we derive equivalent lattice 
gas models at different lattice spacings and determine coefficients directly from binding 
energies and/or few-body scattering data. We also give an estimate of the systematic errors 
and discuss the limits of the lattice model approach in describing long-range ordering. As 
a concrete example we apply the effective field theory lattice gas approach to low-energy 
neutron matter and compare with results from a recent quantum simulation. Despite the 
very low computational cost, essentially the same as that for a 2'^-^-state 3D Ising model, we 
find good agreement with full quantum simulation results when the hopping parameter is 
not too large. 

In our approach temperature-dependent interactions are naturally introduced into the 
lattice gas model. These are necessary to reproduce the physics of the full quantum theory 
at different temperatures. To our knowledge, neutron matter at this density and temperature 
has never been previously described using a lattice gas model. We hope that the broader 
application of lattice gas models as well as the connection with effective field theory will 
prove useful in the study of other many-body fermion systems at nonzero temperature. 

This work was supported in part by the U.S. Department of Energy under grant and 
DE-FG02-04ER41335 and the Minority Graduate Education Program at North Carolina 
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State University. 
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